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Majorana fermions (particles which are their own antiparticle) may or may not exist in Nature 
as elementary building blocks, but in condensed matter they can be constructed out of electron 
and hole excitations. What is needed is a superconductor to hide the charge difference, and a 
topological (Berry) phase to eliminate the energy difference from zero-point motion. A pair of 
widely separated Majorana fermions, bound to magnetic or electrostatic defects, has non-Abclian 
exchange statistics. A qubit encoded in this Majorana pair is expected to have an unusually long 
coherence time. We discuss strategies to detect Majorana fermions in a topological superconduc- 
tor, as well as possible applications in a quantum computer. The status of the experimental search 
is reviewed. 
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I. WHAT ARE THEY? 

A Majorana fermion is a hypothetical particle which is 
its own antiparticle. The search for Majorana fermions 
goes back to the early days of relativistic quantum me- 
chanics. 



A. Their origin in particle physics 

The notion of an antiparticle originated with Paul 
Dirac's 1930 interpretation of the negative-energy solu- 
tions of his relativistic wave equation for spin-| parti- 
cles [1]. The positive-energy solutions describe electrons, 
and the negative-energy solutions correspond to particles 
with the same mass and spin but opposite charge. The 



electron and its antiparticle, the positron, are related by 
a symmetry operation which takes the complex conjugate 
of the wave function. Particle and antiparticle can an- 
nihilate, producing a pair of photons. While the photon 
(described by a real bosonic field) is its own antiparti- 
cle, Dirac fermions are described by complex fields with 
distinct particle and antiparticle. 

In a paper published in 1937, shortly before his disap- 
pearance, Ettore Majorana questioned the need to intro- 
duce a distinct antiparticle for each particle [2] . The com- 
plex Dirac equation can be separated into a pair of real 
wave equations, each of which describes a real fermionic 
field [3] . Majorana suggested that neutral particles might 
be represented by a single real field, and concluded that 
"there is now no need to assume the existence of antineu- 
trons or antincutrinos" . 

We since know that the neutron and antineutron are 
distinct particles, but the neutrino and antineutrino 
could well be the same particle observed in different 
states of motion [4]. It remains to be seen whether or 
not the Majorana fermion will go the way of the mag- 
netic monopole, as a mathematical possibility that is not 
realized by Nature in an elementary particle. 



B. Their emergence in superconductors 

In condensed matter we can build on what Nature of- 
fers, by constructing quasiparticlc excitations with exotic 
properties out of simpler building blocks. This happened 
for magnetic monopoles and it may happen for Majo- 
rana fermions. The strategy to use midgap excitations 
of a chiral p-wave superconductor goes back two decades 
[5-11] (with even earlier traces in the particle physics lit- 
erature [12]). Recent developments in topological states 
of matter have brought this program closer to realization 
[13, 14]. 

The electron and hole excitations of the superconduc- 
tor play the role of particle and antiparticle. Electrons 
(filled states at energy E above the Fermi level) and holes 
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FIG. 1 Profile of the superconducting pair potential A(r) in 
an Abrikosov vortex (solid curve) and bound electron-hole 
states in the vortex core (dashed lines). The left graph shows 
the usual sequence of levels in an s-wave superconductor, ar- 
ranged symmetrically around zero energy. The right graph 
shows the level sequence when superconductivity is induced 
on the surface of a 3D topological insulator, with a nonde- 
generate state at E = 0. This midgap state is a Majorana 
fermion. 
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FIG. 2 Top view of a 2D topological insulator, contacted 
at the edge by two superconducting electrodes separated by 
a magnetic tunnel junction. A pair of Majorana fermions 
is bound by the superconducting and magnetic gaps. The 
tunnel splitting of the bound states depends oc cos(<^>/2) on 
the superconducting phase difference (j>, as indicated in the 
plot. The crossing of the levels at <f> = n is protected by 
quasiparticle parity conservation. 



(empty states at —E below the Fermi level) have opposite 
charge, but the charge difference of 2e can be absorbed as 
a Cooper pair in the superconducting condensate. At the 
Fermi level (E — 0, in the middle of the superconducting 
gap), the eigenstates are charge neutral superpositions of 
electrons and holes. 

That the midgap excitations of a superconductor are 
Majorana fermions follows from electron-hole symmetry: 
The creation and annihilation operators ^(E), ^{E) for 
an excitation at energy E are related by 

1 {E)= 1 \-E), (1) 

At the Fermi level 7(0) = 7 = 7^, so particle and antipar- 
ticle coincide. The anticommutation relation for Majo- 
rana fermion operators has the unusual form 

7«7m + 7m In = 2S nm . (2) 

The operators of two Majoranas anticommute, as for any 
pair of fermions, but the product 7^ = 1 does not vanish. 

Like in the particle physics context, these are just for- 
mal manipulations if the state is degenerate — since a 
Dirac fermion operator a = |(7i + 272) is fully equivalent 
to a pair of Majorana operators. Spin degeneracy, in par- 
ticular, needs to be broken in order to realize an unpaired 
Majorana fermion. The early proposals [5-11] were based 
on an unconventional form of superconductivity, in which 
only a single spin band is involved. Such spin-triplet, 
p-wave pairing is fragile, easily destroyed by disorder. 
Much of the recent excitement followed after Liang Fu 
and Charles Kane showed that conventional spin-singlct, 
s-wave superconductivity could be used, in combination 
with the strong spin-orbit coupling of a topological insu- 
lator [15]. 

The basic mechanism is illustrated in Fig. 1. A three- 
dimensional (3D) topological insulator has an insulating 



bulk and a metallic surface [13, 14]. The 2D surface 
electrons are massless Dirac fermions, very much like in 
graphene — but without the spin and valley degeneracies 
of graphene. A superconductor deposited on the surface 
opens an excitation gap, which can be closed locally by 
a magnetic field. The magnetic field penetrates as an 
Abrikosov vortex, with subgap states E n = (n + a)S, 
n = 0, ±1, ±2, . . ., bound to the vortex core [16]. (The 
level spacing S ~ A 2 /Ep is is determined by the super- 
conducting gap A and the Fermi energy Ep.) Electron- 
hole symmetry restricts a to the values or 1/2. For 
a = the zero-mode Eq = would be a Majorana 
fermion in view of Eq. (1), but one would expect zero- 
point motion to enforce a = 1/2. 

While a = 1/2 indeed holds for the usual massive elec- 
trons and holes, 2D massless Dirac fermions have a = 
— as discovered by Roman Jackiw and Paolo Rossi [12]. 
The reader familiar with graphene may recall the appear- 
ance of a Landau level at zero energy, signifying quanti- 
zation of cyclotron motion without the usual \huj c offset 
from zero-point motion [17]. The absence of a \8 off- 
set in an Abrikosov vortex has the same origin. Massless 
Dirac fermions have their spin pointing in the direction of 
motion. A closed orbit produces a phase shift of tt from 
the 360° rotation of the spin. This Berry phase adds to 
the phase shift of ix in the Bohr-Sommerfeld quantiza- 
tion rule, converting destructive interference at E = 
into constructive interference and shifting the offset a 
from 1/2 to 0. 



C. Their potential for quantum computing 

The idea to store quantum information in Majorana 
fermions originates from Alexei Kitaev [9] . We illustrate 
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the basic idea in Fig. 2 in the context of a 2D topological 
insulator [19, 20], one dimension lower than in Fig. 1. 
The massless Dirac fermions now propagate along a ID 
edge state, again with the spin pointing in the direction of 
motion. (This is the helical edge state responsible for the 
quantum spin Hall effect.) A Majorana fermion appears 
as a zero-mode at the interface between a superconductor 
(S) and a magnetic insulator (I). 

Fig. 2 shows two zero-modes coupled by tunneling in an 
SIS junction, forming a two-level system (a qubit). The 
two states |1) and |0) of the qubit are distinguished by 
the presence or absence of an unpaired quasiparticle. For 
well-separated Majoranas, with an exponentially small 
tunnel splitting, this is a nonlocal encoding of quantum 
information: Each zero-mode by itself contains no infor- 
mation on the quasiparticle parity. 

Dcphasing of the qubit is avoided by hiding the phase 
in much the same way that one would hide the phase 
of a complex number by separately storing the real and 
imaginary parts. The complex Dirac fermion operator 
a = 5(71 + 172) of the qubit is split into two real Ma- 
jorana fermion operators 71 and 72. The quasiparticle 
parity a) a = |(1 + 17172) is only accessible by a joint 
measurement on the two Majoranas. 

While two Majoranas encode one qubit, 2n Majo- 
ranas encode the quantum information of n qubits in 2" 
nearly degenerate states. Without these degeneracies, 
the adiabatic evolution of a state ^ along a closed loop 
in parameter space would simply amount to multiplica- 
tion by a phase factor, h-» e l<x ^>, but now the oper- 
ation may result in multiplication by a unitary matrix, 
<I> 1 y Because matrix multiplications do not com- 

mute, the order of the operations matters. This produces 
the non-Abelian statistics discovered by Gregory Moore 
and Nicholas Read [21], in the context of the fractional 
quantum Hall effect, and by Read and Dmitry Green [8], 
in the context of p-wave superconductors. 

The adiabatic interchange (braiding) of two Majorana 
bound states is a non-Abelian unitary transformation of 
the form 

* 1 Y exp (i^cr z ) *, (3) 

with o~ z a Pauli matrix acting on the qubit formed by the 
two interchanging Majoranas [22, 23]. Two interchanges 
return the Majoranas to their starting position, but the 
final state io~ z ^ is in general not equivalent to the initial 
state <J>. 

An operation of the form (3) is called topological, be- 
cause it is fully determined by the topology of the braid- 
ing; in particular, the coefficient in the exponent is pre- 
cisely 7r/4. This could be useful for a quantum computer, 
even though not all unitary operations can be performed 
by the braiding of Majoranas [24, 25]. 

Before closing this section, we emphasise that the ob- 
ject exhibiting non-Abelian statistics is not the Majorana 
fermion by itself, but the Majorana fermion bound to a 
topological defect (a vortex in Fig. 1, the SI interface 



in Fig. 2, or an e/4 quasiparticle in the fractional quan- 
tum Hall effect). The combined object is referred to as 
an Ising anyon in the literature on topological quantum 
computation [25]. A free Majorana fermion (such as may 
be discovered in particle physics) has ordinary fermionic 
statistics — it is not an Ising anyon. The same applies 
to unbounded Majorana fermions at the edge or on the 
surface of a topological superconductor. 

In what follows we will concentrate on the Majoranas 
bound to a topological defect, because of their exotic 
statistics. We could have included the Ising anyons in the 
fractional quantum Hall effect, but because that topic is 
already very well reviewed [26] , we limit ourselves to the 
superconducting implementations. 



II. HOW TO MAKE THEM 

The route to Majorana fermions in superconductors 
can follow a great variety of pathways. The growing list 
of proposals includes Refs. 5-11, 15, 18-20, 27-60. There 
are so many ways to make Majorana fermions because 
the requirements are so generic: Take a superconduc- 
tor, remove degeneracies by breaking spin-rotation and 
time-reversal symmetries, and then close and reopen the 
excitation gap. As the gap goes through zero, Majo- 
rana fermions emerge as zero-modes bound to magnetic 
or electrostatic defects [61, 62]. We summarize the main 
pathways, and refer to a recent review [63] for a more 
detailed discussion. 



A. Shockley mechanism 

From this general perspective, Majorana bound states 
can be understood as superconducting counterparts of 
the Shockley states from surface physics [64, 65]. The 
closing and reopening of a band gap in a chain of atoms 
leaves behind a pair of states in the gap, bound to the 
end points of the chain, see Fig. 3. Shockley states are 
unprotected and can be pushed out of the band gap by 
local perturbations. In a superconductor, in contrast, 
particle-hole symmetry requires the spectrum to be ±E 
symmetric, so an isolated bound state is constrained to 
lie at E = and cannot be removed by any local pertur- 
bation, see Fig. 4. 

The closing of the excitation gap, followed by its re- 
opening with opposite sign, is a topological phase tran- 
sition. The phase transition is called topological, in dis- 
tinction to thermodynamic, because the sign Q = ±1 
of the gap cannot be seen in thermodynamic properties. 
This so-called Z 2 topological quantum number counts the 
parity of the number of Majorana fermions bound to the 
defect. Only an odd number of Majoranas (Q = — 1) 
produces a stable zero-mode. A defect with Q = —1 is 
called topologically nontrivial, while for Q = +1 it is 
called topologically trivial. 

Because the Majorana fermions are constructed from 
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FIG. 3 Illustration of the Shockley mechanism for the for- 
mation of bound states at the end points of an atomic chain. 
The lower panel shows the potential profile along the chain 
and the upper panel shows the corresponding energy levels 
as a function of the atomic separation a. The end states ap- 
pear upon the closing and reopening of the band gap. Figure 
adapted from Ref. 64. 




FIG. 4 Emergence of a pair of zero-energy Majorana states 
in a model calculation of a chiral p-wave superconductor con- 
taining an electrostatic line defect. The gap closes and re- 
opens as the defect potential Uo + SU is made more and more 
negative, at fixed positive background potential Uo- The in- 
set shows the probability density of the zero-mode in the 2D 
plane of the superconductor, with the line defect along the 
y-axis. Figure adapted from Ref. 65. 



ordinary Dirac electrons, an unpaired Majorana at a 
topologically nontrivial defect must have a counterpart 
somewhere else in the system. The two Majoranas are 
evident in Figs. 2 and 4. In Fig. 1 the second Majorana 
extends along the outer perimeter of the superconduc- 
tor. One could try to eliminate this second Majorana 
by covering the entire topological insulator by a super- 
conductor. But then the flux line would intersect the 
superconductor at two points, producing again a pair of 



FIG. 5 The red solid curves locate the closing of the excita- 
tion gap of an electrostatic line defect (width W = 4ft/mA) 
in a chiral p-wave superconductor, described by the Hamil- 
tonian (4). In the shaded regions the defect is topologically 
nontrivial (Q — —1), with Majorana states bound to the two 
ends. Figure adapted from Ref. 65. 

Majoranas. 

Let us see how these topological phase transitions ap- 
pear in some representative systems. 

B. Chiral p-wave superconductors 

The closing and reopening of the gap in Fig. 4 is de- 
scribed by the Bogoliubov-De Gennes Hamiltonian of a 
2D chiral p-wave superconductor, 

H _ ( U + p 2 /2m A(p x - i Py ) \ , . 

\A*(p x +i Py ) -U-p 2 /2mJ- w 

The diagonal elements give the electrostatic energy 
±U(r) and kinetic energy ±p 2 /2m of electrons and holes. 
(Energies are measured relative to the Fermi level.) The 
off-diagonal elements couple electrons and holes via the 
superconducting pair potential, which has the chiral p- 
wave orbital symmetry oc p x ±ip y . (For equal-spin triplet 
pairing the spin degree of freedom can be omitted.) 

The line defect is constructed by changing U from the 
background value Uq to Uo + SU in a strip of width W. 
As SU is varied, multiple closings and reopenings of the 
gap appear, see Fig. 5. This is the 2D generalization 
[65-68] of the ID Kitaev chain [9]. The gap closing is a 
result of destructive interference of transverse modes in 
the strip. Each new mode is associated with one closing- 
reopening of the gap, so that the defect is topologically 
nontrivial (Q = —1) for an odd number of modes and 
trivial (Q = +1) for an even number of modes. 

Electrostatic line defects are one way of producing Ma- 
jorana fermions in a chiral p-wave superconductor. Mag- 
netic vortices are another way [6, 8, 69-72], those defects 
are topologically nontrivial for U < 0. Strontium ruthen- 
ate (Sr 2 Ru04) is a candidate p-wave material to observe 
the predicted zero modes [11]. 
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FIG. 6 Dirac-to-Majorana fermion converter on the surface 
of a 3D topological insulator. Arrows indicate the propat- 
ing modes at the interface between a superconductor and a 
magnetic insulator and at the magnet-magnet interface. An 
electron (Dirac fermion) injected by the voltage source at the 
right is split into a pair of Majorana fermions. These fuse at 
the left, either back into an electron or into a hole, depend- 
ing on whether the superconductor contains an even or an 
odd number of vortices. The recombination as a hole adds a 
Cooper pair to the superconductor, which can be detected in 
the current to ground. Figure adapted from Ref. [74]. 

C. Topological insulators 

In a topological insulator the closing and reopening of 
the band gap is a consequence of strong spin-orbit cou- 
pling, which inverts the order of conduction and valence 
bands [13, 14]. The surface of a topological insulator 
supports nondegenerate, massless Dirac fermions, with 
Hamiltonian 

H = v F p ■ a + U + Ma z . (5) 

A 3D topological insulator, such as Bi2Se3 or Bi2Te3, has 
2D Dirac fermions on the surface, while a 2D topological 
insulator, such as a HgTe/CdTe or InAs/GaSb quantum 
well, has ID Dirac fermions along the edge. The term per 
represents p x a x +p v <J y or p x a x for surface or edge Dirac 
fermions, respectively. The extra term Ma z accounts for 
the exchange energy from a magnetic insulator. 

As illustrated in Figs. 1 and 2, both the surface and 
edge Dirac fermions can give rise to Majorana bound 
states [15]. The Bogoliubov-De Gennes Hamiltonian that 
describes these zero modes has the form 

W£° % )■ (6) 

The diagonal contains the electron and hole Dirac Hamil- 
tonians (5), related to each other by the time-reversal 
operation H t— > ayH^ay (which inverts p and er). The 
pair potential A (which can be complex as a result of a 
magnetic vortex) is induced by s-wave superconductivity, 
so it is momentum independent. 

Fig. 6 illustrates the conversion of Dirac fermions into 
Majorana fermions on the surface of a 3D topological 



insulator [73, 74]. Part of the surface is covered with 
a superconductor (M = 0,A / 0) and part is covered 
with a magnetic insulator (A = 0, M ^ 0). Both A 
and M open a gap for the surface states. The gap closes 
at the M-A interface and also at the M^-M^ interface 
between opposite magnetic polarizations (from M-f > 
to Mi < 0). The ID interface states propagate only in 
a single direction, much like the chiral edge states of the 
quantum Hall effect. 

The Mf—M± interface leaves electrons and holes un- 
coupled, so there are two modes at the Fermi level, one 
containing electrons and one containing holes. These are 
Dirac fermions, with distinct creation and annihilation 
operators a* and a. The electron-hole degeneracy is bro- 
ken at the M-A interface, which supports only a sin- 
gle mode 7 = 7^ at the Fermi level. The M^-M^-A 
tri-j unction splits an electron or hole into two Majorana 
fermions, 

|(7i + *72)> a f -> §(71 - ry 2 ). (7) 

The inverse process, the fusion of two Majorana 
fermions into an electron or hole, happens after the Ma- 
jorana fermions have encirled the superconductor and 
picked up a relative phase shift. At the Fermi level this 
phase shift is entirely determined by the parity of the 
number n of vortices in the superconductor. For even n 
the fusion conserves the charge, while for odd n a charge 
2e is added as a Cooper pair to the superconductor (so an 
electron is converted into a hole and vice versa) . The two 
processes can be distinguished experimentally by measur- 
ing the current to ground in the superconductor. Since 
each vortex binds one Majorana fermion (see Fig. 1), the 
structure of Fig. 6 can be seen as an electrical interfer- 
ometer in which mobile Majoranas measure the parity of 
the number of Majorana bound states. 



D. Semiconductor heterostructures 

The gap inversion in a topological insulator happens 
without superconductivity. One might alternatively try 
to directly invert the superconducting gap. In an s-wave 
superconductor a magnetic field closes the gap, but how 
to reopen it? A promising strategy is to rely on the 
competing effects of a spin-polarizing Zeeman energy and 
a depolarizing spin-orbit coupling [29, 30, 75]. 

The 2D electron gas of a semiconductor heterostruc- 
ture, such as an InAs quantum well, has a strong spin- 
orbit coupling from the Rashba effect. The orbital effect 
of a parallel magnetic field is suppressed, leaving only the 
spin-polarizing Zeeman effect. These two effects compete 
in the Hamiltonian 

2 

H ° = 2m +U + ^ir( axp y ~ a y px ^ + \9cS^bBg x . (8) 

Characteristic length and energy scales are l BO = 
h 2 /ma so and E so = ma 2 /h 2 . Typical values in InAs 
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FIG. 7 Closing and reopening of the band gap in the Hamiltonian (8) of a 2D semiconducting nanowire (width W = l so ) on 
a superconducting substrate in a parallel magnetic field (A = 10-E S o, Ez = 10.5 E BO ). The eight panels show the excitation 
energy near the Fermi level (E — 0) as a function of the wave vector k along the nanowire, for different values of the chemical 
potential n = —U (listed in units of the spin-orbit coupling energ E so ). The colors blue or yellow of the panels indicate that the 
system is in a topologically trivial or nontrivial phase, respectively. The topological phase transition occurs in the uncolored 
panels. The nanowire supports Majorana bound states in the yellow panels. Data supplied by M. Wimmer. 



are l so — 100 nm, E so = 0.1 meV. The Zeeman en- 
ergy is Ez — \g c sl^BB = lmeV at a magnetic field 
B = 1 T. A superconducting proximity effect with a 
type-II superconductor like Nb is quite possible at these 
field strengths. The pair potential A induced in the 2D 
electron gas then couples electrons and holes via the 
Bogoliubov-De Gennes Hamiltonian (6) [now with Hq 
given by Eq. (8)]. 

As discovered in Refs. 32, 33, the resulting band gap in 
a nanowire geometry closes and reopens upon variation 
of electron density (through a variation of U) or magnetic 
field (see Fig. 7). Majorana bound states at the two ends 
of the nanowire alternatingly appear and disappear at 
each of these topological phase transitions. 



III. HOW TO DETECT THEM 

Majorana fermions modify the transport properties 
and the thermodynamic properties of the superconduc- 
tor, providing ways to detect them. We summarize some 
of the signature effects of Majorana fermions, others can 
be found in Refs. 20, 73, 74, 76-87. 




Ef - Vqpc [E so ] 

FIG. 8 Solid curves: conductance of a ballistic NS junction, 
with the superconductor in a topologically trivial or non- 
trivial phase. The dotted curve is for an entirely normal sys- 
tem. The data is calculated from the model Hamiltonian (8) . 
The point contact width is varied by varying the potential 
Vqpc inside the constriction at constant Fermi energy Ef- 
The dotted horizontal lines indicate the shift from integer to 
half-integer conductance plateaus upon transition from the 
topologically trivial to nontrivial phase. Figure adapted from 
Ref. 92. 



A. Half-integer conductance quantization 

Tunneling spectroscopy is a direct method of detec- 
tion of a Majorana bound state [88, 89]: Resonant tun- 
neling into the midgap state produces a conductance of 
2e 2 /h, while without this state the conductance vanishes. 
A complication in the interpretation of tunneling spec- 



troscopy is that the zero-bias peak may be obscured by 
resonances from subgap states at nonzero energy [90, 91]. 
A ballistic point contact provides a more distinctive sig- 
nature of the topologically nontrivial phase [92] , through 
the half-integer conductance plateaus shown in Fig. 8. 

Both the tunneling and ballistic conductances can be 
understood from the general relation [93] between the 
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conductance G of a normal-metal-superconductor (NS) 
junction and the Andreev reflection eigenvalues i?„, 



G 



h 



22 Rr, 



(9) 



The RnS represent the probability for Andreev reflection 
in the n-th eigenmode at the Fermi level. The factor of 
two is not due to spin (which is included in the sum 
over n) , but due to the fact that Andreev reflection of an 
electron into a hole doubles the current. 

There is no time-reversal symmetry, so Kramers de- 
generacy does not apply. Still, particle-hole symmetry 
requires that any R n is twofold degenerate (Beri de- 
generacy [94]) — with two exceptions: R n = and 
R n = f may be nondegenerate. The nondegenerate 
Andreev reflection eigenvalue from a Majorana bound 
state is pinned to unity, contributing to the conductance 
a quantized amount of 2e 2 jh. All other fully Andreev 
reflected modes are twofold degenerate and contribute 
4e 2 //i. The resulting conductance plateaus therefore ap- 
pear at integer or half-integer multiples of 4e 2 //i, depend- 
ing on whether the superconductor is topologically trivial 
or not. 

The plateaus at (n + 1/2) x 4e 2 //i are reminiscent of 
the quantum Hall plateaus in graphene, and both orig- 
inate from a zero mode, but the sensitivity to disorder 
is entirely different. The topological quantum number 
Q £ Z for the quantum Hall effect, while Q £ Z 2 for a 
topological superconductor. The corresponding topolog- 
ical protection against disorder extends to all plateaus 
for the quantum Hall effect, but only to the lowest n = 
plateau for the topological superconductor. 

It may appear paradoxical [95, 96] to have an elec- 
trical current flowing through a single Majorana bound 
state, since one Majorana fermion operator 7 represents 
only half of an electronic state. However, the Hermitian 
operator 2(0 + 0^)7 is a local coupling of Dirac and Majo- 
rana operators [88] , so electrical conduction can be a fully 
local process — involving only one of the two spatially 
separated Majorana fermions. 
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FIG. 9 Majorana bound states (red) at the edge of a 2D topo- 
logical insulator (cf. Fig. 2), split into a pair of levels at ±Em 
by a nonzero overlap. The levels are broadened due to a tun- 
nel coupling ri,F2 through the magnet to the outside edge 
state. An electron incident from the left on the grounded su- 
perconductor can be Andreev reflected as a hole, either locally 
(to the left) or nonlocally (to the right). Nonlocal Andreev 
reflection is equivalent to the splitting of a Cooper pair by the 
two Majoranas. For Fi,r2 <C Em local Andreev reflection is 
suppressed. Figure adapted from Ref. 20. 



presence or absence of an unpaired quasiparticle, see 
Section I.C.) 

Nonlocal tunneling processes appear if the level split- 
ting Em is large compared to the level broadening Ti, T 2 . 
For a grounded superconductor the nonlocality takes the 
form of nonlocal Andreev reflection, which amounts to a 
splitting of a Cooper pair by the two Majorana bound 
states [20] (see Fig. 9). The Cooper pair splitting can 
be detected in a noise measurement through a positive 
cross-correlation of the currents I\ and I 2 to the left and 
right of the superconductor. 

For an electrically isolated superconducting island any 
charge transfer onto the island is forbidden by the charg- 
ing energy, so there can be no Andreev reflection. An 
electron incident on one side of the island is either re- 
flected to the same side or transmitted, still as an elec- 
tron, to the other side. The nonlocality [97] now appears 
in the ratio of the reflection and transmission probabili- 
ties on resonance, 



r/t = {v 1 -r 2 ) 2 /(r 1 + r 2 ) 2 



(10) 



B. Nonlocal tunneling 

Nonlocal conduction involving both Majorana 
fermions becomes possible if there is a coupling between 
them [20, 89, 97-102]. The coupling term has the 
generic form i-Ej\/7i72i with eigenvalues ±Em- The 
energy Em may be a tunnel coupling due to overlap of 
wave functions, in which case it decays exponentially 
oc e~ d /^° with the ratio of the separation d of the 
Majoranas and the superconducting coherence length 
£o- If the superconductor is electrically isolated (not 
grounded) and of small capacitance C, then the charging 
energy Em — e 2 /C provides a Coulomb coupling even 
without overlap of wave functions. (Recall that the two 
states of a pair of Majoranas are distinguished by the 



which is independent of the size of the island. No matter 
how far the two Majoranas are separated, the charging 
energy couples them into a single electronic level. In par- 
ticular, for identical tunnel couplings Ti = r 2 the elec- 
tron is resonantly transmitted through the island with 
unit probability . 



C. 47r-periodic Josephson effect 

So far we discussed signatures of Majoranas in the elec- 
trical conduction out of equilibrium, in response to a volt- 
age difference between the superconductor and a normal- 
metal electrode. In equilibrium an electrical current (su- 
percurrent) can flow between two superconductors in the 
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absence of any applied voltage. This familiar DC Joseph- 
son effect [103] acquires a new twist [9, 19, 32, 33, 104- 
113] if the junction between the superconductors contains 
Majorana fermions, as in Fig. 2. 

Quite generally, the supercurrent Ij is given by the 
derivative 



of the energy E of the Josephson junction with respect 
to the superconducting phase difference <f>. While in the 
conventional Josephson effect only Cooper pairs can tun- 
nel (with probability t<1), Majorana fermions enable 
the tunneling of single electrons (with a larger probabil- 
ity y/r). The switch from 2e to e as the unit of trans- 
ferred charge between the superconductors amounts to a 
doubling of the fundamental periodicity of the Josephson 
energy, from E oc cos</> to E oc cos(0/2). 

If the superconductors form a ring, enclosing a flux $, 
the period of the flux dependence of the supercurrent Ij 
doubles from 2ir to 4-7T as a function of the Aharonov- 
Bohm phase 2eQ/h. This is the 47r-periodic Josephson 
effect [9, 104]. As a function of the enclosed flux, Ij 
has the same h/e periodicity as the persistent current 
In through a normal-metal ring (radius L), but the size 
dependence is entirely different: While In decays asl/L 
or faster, Ij has the //-independence of a supercurrent. 

Since the two branches of the E-<j> relation differ by 
one unpaired quasiparticle (see Fig. 2), external tun- 
neling events which change the quasiparticle parity (so- 
called quasiparticle poisoning) restore the conventional 
27r-periodicity [19] . In a closed system, the 47r-periodicity 
is thermodynamically stable, provided that the entire 
ring is in a topologically nontrivial state (to prevent 
quantum phase slips) [109]. 

D. Thermal metal-insulator transition 

Collective properties of Majorana fermions can be de- 
tected in the thermal conductance. Superconductors are 
thermal insulators, because the excitation gap A sup- 
presses the energy transport by quasiparticle excitations 
at low temperatures To -C A/ks- Disorder can create 
states in the gap, but these are typically localized. How- 
ever, the Majorana midgap states in a topological super- 
conductor can give rise to extended states, since they are 
all resonant at the Fermi level. This transforms a thermal 
insulator into a thermal metal [114, 115]. 

The thermal metal-insulator transition is called a class 
D Anderson transition, in reference to a classification of 
disordered systems in terms of the presence or absence 
of time-reversal, spin-rotation, and particle-hole symme- 
try [116-118]. Class D has only particle-hole symmetry. 
The chiral p-wave superconductor is a two-dimensional 
system in class D. Its thermal transport properties in the 
absence of Majorana fermions are similar to the electri- 
cal transport properties in the quantum Hall effect (class 
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FIG. 10 Main plot and right inset: Average density of states 
p in a model calculation of a chiral p-wave superconductor 
with electrostatic disorder [65]. The Hamiltonian (4) is dis- 
cretized on a lattice of size 400 a x 400 a and the potential 
fluctuates randomly from site to site (r.m.s. AU). Majorana 
fermions produce a midgap peak in the density of states. Left 
inset: Average thermal conductivity a in a strip geometry of 
length L and width W = 5L, for the same Hamiltonian (4) 
but calculated with a different method of discretization [119]. 
Data points of different color correspond to different disor- 
der strengths AU and different scattering lengths £. Upon 
increasing disorder, a transition from insulating to metallic 
scaling is observed. In the metallic phase, the conductivity 
and density of states have a logarithmic dependence on, re- 
spectively, system size and energy. 



A, all symmetries broken). The bulk is insulating while 
the boundary supports chiral (unidirectional) edge states 
that give rise to the thermal quantum Hall effect [7, 8]. 
The thermal analogue of the conductance quantum e 2 /h 
is Go = T: 2 k 2 B T a /6h. 

The correspondence between thermal and electrical 
quantum Hall effect breaks down in the presence of Majo- 
rana fermions. Their collective effect is illustrated in Fig. 
10, obtained from the Hamiltonian (4) of a chiral p-wave 
superconductor [65, 119]. (Similar results have been ob- 
tained in other models of Majorana fermions [120-122].) 
A randomly varying electrostatic potential creates a ran- 
dom arrangement of Majorana midgap states, via the 
Shockley mechanism of Fig. 4. The states are slightly 
displaced from E — by the overlap of wave functions. 
The resulting density of states has the logarithmic profile 
cx In \E\, responsible for the logarithmic size dependence 
of the thermal conductivity [7] , 



a = (Gq/tt) InL + constant. (12) 



The thermal metal-insulator transition of Fig. 10 has no 
electrical analogue. 
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IV. HOW TO USE THEM 

Finding Majorana fermions in a superconductor is re- 
warding in and of itself. These particles might also 
provide a fundamentally new way to store and manip- 
ulate quantum information, with possible applications in 
a quantum computer. 

A. Topological qubits 

In Section I.C we considered a qubit formed out of a 
pair of Majorana fermions. The two states |0) and |1) of 
this elementary qubit differ by quasiparticle parity, which 
prevents the creation of a coherent superposition. For a 
quantum computation we combine two elementary qubits 
into a single logical qubit, consisting of four Majorana 
fermions [25] . Without loss of generality one can assume 
that the joint quasiparticle parity is even. The two states 
of the logical qubit are then encoded as 1 00) and |11). 
These two states have the same quasiparticle parity, so 
coherent superpositions are allowed. 

An arbitrary state |'J / ) of the logical qubit has the form 

|*)=a|00)+ J 8|ll) 1 |a| 2 + |/3| 2 = l. (13) 

Pauli matrices in the computational basis |00), |11) are 
bilinear combinations of the four Majorana operators, 

<7 X = -«7273 o-y = 17173, <J Z = -17172. (14) 

It is said that the qubit (13) is topologically protected 
from decoherence by the environment [9], because the 
bit-flip or phase-shift errors produced by the Pauli ma- 
trices (14) can only appear if there is a coupling between 
pairs of Majorana fermions. The two types of coupling 
were discussed in Section III.B: tunnel coupling when 
the Majoranas are separated by less than a coherence 
length, and Coulomb coupling when the Majoranas are 
on a superconducting island of small capacitance. 

The topological protection relies on the presence of a 
nonzero gap for quasiparticle excitations. Sub-gap exci- 
tations may exchange a quasiparticle with a Majorana 
fermion, provoking a bit-flip error. Error correction is 
possible if the sub-gap excitations remain bound to the 
Majorana fermion, in particular sub-gap excitations in 
a vortex core are not a source of decoherence [123-125]. 
The topological protection does not apply if the super- 
conductor is contacted by a gapless metal, allowing for 
the exchange of unpaired electrons (quasiparticle poison- 
ing) [126, 127]. 

B. Read out 

To read out a topological qubit one needs to remove the 
topological protection by coupling the Majorana fermions 
and then measure the quasiparticle parity. Tunnel cou- 
pling is one option, for example in the geometry of Fig. 2 



Cooper pair box 




topological qubit 



FIG. 11 Read out of a topological qubit in a Cooper pair 
box. Two superconducting islands (red), connected by a 
split Josephson junction (crosses) form the Cooper pair box. 
The topological qubit is formed by two pairs of Majorana 
fermions (yellow dots), at the end points of two undepleted 
segments (blue) of a semiconductor nanowire (shaded ribbon 
indicates the depleted region) . A magnetic fiux $ enclosed by 
the Josephson junction controls the charge sensitivity of the 
Cooper pair box. To read out the topological qubit, one pair 
of Majorana fermions is moved onto the other island. Depend- 
ing on the quasiparticle parity, the resonance frequency in a 
superconducting transmission line enclosing the Cooper pair 
box (green) is shifted upwards or downwards by the amount 
given in Eq. (16). Figure adapted from Ref. 130. 

the quasiparticle parity of two Majorana fermions can be 
measured by the difference in tunnel splitting [19]. The 
alternative Coulomb coupling allows the joint read out of 
more than a single qubit [128-130]. We concentrate on 
this second option, since two-qubit read out is required 
for quantum computations [25]. 

Consider a superconducting island (charge Q, super- 
conducting phase 0), containing 2N Majorana fermions. 
Cooper pairs can enter and leave the island via a Joseph- 
son junction. The read-out operation amounts to a mea- 
surement of Q modulo 2e. The conventional even-odd 
parity effect of the superconducting ground state does 
not apply here, because there is no energy cost of A 
for an unpaired electron in a midgap state. Indeed, the 
quasiparticle parity V does not enter in the Hamiltonian 
directly, but as a constraint on the eigenstates [97, 131], 

= e i7r7 >(<£ + 2tt), V= i + ^ A '7i72---72iv. (15) 

This constraint enforces that the eigenvalues of the 
charge operator Q = —2ei d/d<j> are even or odd multiples 
of e when V equals or 1, respectively. 

The parity constraint (15) modifies the energy of 4> > 
(f>+2jT quantum phase slips of the superconducting island, 
induced by the nonzero charging energy Eq — e 2 /2C. 
The V dependence can be measured spectroscopically in 
a SQUID geometry [128] or in a Cooper pair box [130]. 
We show the latter geometry in Fig. 11. 

The magnitude of the P-dependent energy shift AE 
in the Cooper pair box is exponentially sensitive to the 
ratio of Josephson and charging energies, 

AE = (2P-1)U, U~ y/huj p Eje- R ^ /Ej . (16) 

(The frequency uj p — ^/8EqE j Jh is the Josephson 
plasma frequency.) By varying the flux $ through a 
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split Josephson junction, the Josephson energy Ej oc 
cos(e<I>//i) becomes tunable. In the transmon design of 
the Yale group, a variation of Ej/Eq over two orders 
of magnitude has been realized [132]. The Coulomb 
coupling U of the Majorana fermions can therefore be 
switched on and off by varying the flux. 

C. Braiding 

In the two-dimensional geometry of Fig. 1 the Majo- 
rana bound states can be exchanged by moving the vor- 
tices around [15]. The Majorana fermions in the one- 
dimensional geometry of Fig. 11 are separated by insu- 
lating regions on a single nanowire, so they cannot be 
exchanged (at least not without rotating the wire it- 
self [133]). The exchange of Majorana fermions, called 
"braiding", is needed to demonstrate their non-Abelian 
statistics [8] . It is also an essential ingredient of a topo- 
logically protected quantum computation [24]. In order 
to be able to exchange the Majoranas one can use a sec- 
ond nanowire, running parallel to the first and connected 
to it by side branches [134, 135]. 

The minimal Hamiltonian that can describe the braid- 
ing contains three Majorana fermions 71,72,73 coupled 
to a fourth one 70, 

3 

H = J2 U kijolk- (17) 
fc=i 

The three parameters Uf. > can describe tunnel cou- 
pling [136] (tunable by a gate voltage) or Coulomb cou- 
pling [137] (tunable by the flux through a Josephson junc- 
tion). A tri-junction of three Cooper pair boxes that is 
described by this Hamiltonian is shown in Fig. 12. 

The braiding operation consists of three steps, de- 
noted C>3i, 0i2, and 023- At the beginning and at 
the end of each step two of the couplings are off and 
one coupling is on. The step O^k' consists of the se- 
quence { k, k ' } = { on, off} i-» { on, on} h->- { off, on} . The ef- 
fect of this sequence is to transfer the uncoupled Majo- 
rana 7^' 1 — ^ — 7fe. (The minus sign appears in order to 
conserve the quasiparticle parity.) The result after the 
three steps shown in Fig. 12 is that the Majoranas at sites 
1 and 2 are switched, with a difference in sign, 72 <— > —71, 
7i ^72- The corresponding adiabatic time evolution op- 
erator in the Heisenberg representation 7^ i— > U^ftfA is 
given by 

U = (I + 7172) = cxp (^7172) = exp (ij°-z) ■ (18) 

This is the operator of Eq. (3), representing a non- 
Abelian exchange operation. 

V. OUTLOOK ON THE EXPERIMENTAL PROGRESS 

As we have seen in Section II, there is no shortage 
of proposals for superconducting structures that should 
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FIG. 12 Lower panel: Three Cooper pair boxes connected at 
a tri-junction via three overlapping Majorana fermions (which 
effectively produce a single zero-mode 70 = 3 _1//2 (7i+7 2 + 73) 
at the center). Upper panel: Schematic of the three steps 
of the braiding operation. The four Majoranas of the tri- 
junction (the three outer Majoranas 71,72,73 and the effec- 
tive central Majorana 70) are represented by circles and the 
coupling Uk is represented by lines (solid in the on state, 
dashed in the off state). White circles indicate strongly cou- 
pled Majoranas, colored circles those with a vanishingly small 
coupling. The small diagram above each arrow shows an in- 
termediate stage, with one Majorana delocalized over three 
coupled sites. The three steps together exchange the Majo- 
ranas 1 and 2, which is a non-Abelian braiding operation. 
Figure adapted from Ref. 137. 



bind a Majorana zero-mode to a magnetic vortex or elec- 
trostatic defect. This gives much hope for a variety of 
experimental demonstrations in the coming years. There 
has already been a remarkable progress. 

A Josephson effect at the surface of a 3D topologi- 
cal insulator with superconducting electrodes has been 
observed in BiSb alloys [138], and in crystalline Bi 2 Sc 3 
[139-142] and Bi 2 Te 3 [143, 144]. These experiments, and 
related Andreev conductance measurements [145-148], all 
involve wide electrodes with a macroscopic number of oc- 
cupied modes at the Fermi level. While the Josephson 
effect and the Andreev conductance show interesting and 
unusual features, these cannot be readily attributed to 
the single Majorana zero-mode (typically only one out 
of 10 5 modes). Vortices (as in Fig. 1), or other means 
of confinement would be needed to produce a Majorana 
bound state. 

The edge of a 2D topological insulator provides a 
single-mode conductor that could support spatially sep- 
arated Majorana bound states, as in Fig. 2. A su- 
perconducting proximity effect has been observed in an 
InAs/GaSb quantum well [149], and also HgTe/CdTe 
would be a promising system — if the Majorana can be 
confined to the superconducting interface by a magnetic 
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FIG. 13 Top photograph: InSb wire between a normal-metal 
(N) and a superconducting (S) electrode. A barrier gate cre- 
ates a confined region (marked in red) at the interface with 
the superconductor. Other gates are used to locally vary the 
electron density. A magnetic field B is applied parallel to the 
wire. Bottom graph: Differential conductance at 60 mK for B 
incrementing from to 490 mT in 10 mT steps. (Traces are 
offset for clarity, except for the lowest trace at B — 0.) The 
peaks at ±250 /xeV correspond to the gap induced in the wire 
by the superconducting proximity effect. Upon increasing B 
a peak develops at zero voltage, signaling the appearance of a 
Majorana zero-mode in the confined region. Figure adapted 
from Ref. 150. 



insulator or magnetic field. 

At this time of writing (April 2012), semiconductor 
nanowires have come furthest in the realization of Ma- 
jorana fermions, following the proposal of Lutchyn et al. 
[32] and Oreg et al. [33]. Convincing evidence for a Ma- 
jorana zero-mode in an InSb nanowire has been reported 
by Kouwenhoven and his group [150], see Fig. 13. These 
developments give hope that the rich variety of unusual 
properties of Majorana fermions, reviewed in this article, 
will soon be observed experimentally. 
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